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We present a novel and convenient mean-field method, and apply it to study the metal- 
lic/antiferromagnetic interface of ft:-(BEDT-TTF)2Cu[N(CN)]2Y organic superconductors (BEDT- 
TTF is bis-ethylendithio-tetrathiafulvalene, Y=Cl,Br). The method, which fully exploits the crystal 
symmetry, allows one to obtain the mean-field solution of the two-dimensional Hubbard model for 
very large lattices (tipically 6x10^ sites), yielding a reliable description of the phase boundary in 
a wide region of the parameter space. The metal/antiferromagnet transition appears to be second 
order, except for a narrow region of the parameter space, where the transition is very sharp and 
possibly first order. The coexistence of metallic and antiferromagnetic properties is only observed 
for the transient state in the case of smooth second order transitions. The relevance of the present 
results to the complex experimental behavior of centrosymmetric k-(BEDT-TTF)2Cu[N(CN)]2Y 
salts is discussed. 
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I. INTRODUCTION 

The K-phase (BEDT-TTF)2X sahs exhibit a great va- 
riety of physical properties as a function of tempera- 
ture, pressure, anion (X) substitution, deuteration, and 
even disorder in the ethylene end-groups. Supercon- 
ducting (SC), antiferromagneti&i (AF), metallic and in- 
sulating phases are observed.El'a Of particular interest 
is the AF/SC/metal borderhne, which for k-(BEDT- 
TTF)2Cu(NCS)2 and k-(BEDT-TTF)2Cu[N(CN)]2Y (Y 
= CI, Br; hereafter ET-Y family) occurs in a very narrow 
region of the temperature-pressure (T,p) space. For the 
aforementioned compounds, a schematic zero tempjer, 
ture phase diagram can be drawn as shown in fig. 
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FIG. 1. Universal zero temperature phase diagram for 
K-phase BEDT-TTF salts. 



The proximity of SC and AF phases, together with 
other experimental evidences, suggested a possible roli 
spin fluctuations in the superconductivity mechanism 
and prompted inteaspe theoretical investigation on the 
SC/AF borderline.0ill|-Since the early suggestion by 
Kino and FukuyamaJ3't3 mean-field (mf) approaches 
have often, Jaeen adopted to investigate the SC/AF 
interface.H'u'tj Several previous,jp£ treatments consid- 
ered low-symmetry structures, fflllj namely k-(BEDT- 
TTF)2Cu(NCS)2, with four inequivalent molecules in 
the unit cell. The resulting numerical calculation is 
complex and computationally very demanding, so that 
only fairly small lattices have been considered, leading 
to large intrinsic uncertainties on the estimated proper- 
ties, particularly at the phase transition. On the other 
hand, in the orthorhombic centrosymmetric K-phase crys- 
tals of thCj-ET-Y family all molecules in the layer are 
equivalent .E3o The greater symmetry with respect to 
K-(BEDT-TTF)2Cu(NCS)2 apparently does not lead to 
significant differences in the physical behavior, notably 
in the SC properties. Indeed, ET-Br is a supercon- 
ductor at ambient pressure, and ET-Cl under moder- 
ate pressure presents the highest Tc observed in (BEDT=. 
TTF)2X salts.El The mf approach has been applieda'Eil 
also to these more symmetric lattices to study the co- 
existence of charge ordering instabilities. Including the 
order parameters relevant to charge instabilities further 
increases the complexity of the calculation and leads 
to large uncertainties in the transition region, as rec- 
ognized by the authors of Ref. ^. We instead propose 
to fully exploit the higher symmetry of the centrosym- 
metric ET-Y structures to build a symmetrized mf ap- 



proach for the supercoiidpcting metallic/magnetically or- 
dered phase transition.E3 By exploiting symmetry we are 
able to work on lattices as large as 6x10^ sites, and to 
explore a wide region of the parameter space, keeping 
numerical procedures and finite-size effects under con- 
trol. Moreover and most importantly, we get a simple 
and complete description of the electronic bands of these 
system, focusing only on those effects which are directly 
connected with the interesting physics. The mf treat- 
ment maps the problem of interacting electrons into an 
effective non-interacting Hamiltonian. The comparison 
between reliable mf results and available experimental 
data then allows us to safely define the intrinsic limita- 
tions of effective one-electron pictures in describing the 
physics of K-phase salts. 

Recently more refined approaches have been ap- 
plied to K-phase a a lts. For instance, fluctuation- 
exchange^FLEX)E3t3 or third-order perturbation (PT) 
approachtZl have been applied to investigate the fluctu- 
ation mechanism for superconductivity and to estimate 
the critical temperature. However these approaches are 
only valid in the weak electronic correlations regime, and 
have been applied to a simpliffed-jpcidol for the K-phase 
layer, the so-called dimer model.Bl3^ESE3 The same model 
has. also been adopted in the framework of the dynamical 
mfl'clj and within a renormalization group approachjEll 
which offer complementary information with respect to 
ordinary mf techniques. As we shall discuss in the follow- 
ing, the reliability of the dimer model cannot be taken 
for granted in the whole parameter space. Here we show 
that by exploting the high symmetry of the ET-Y salts 
family one can get a picture of the K-phase layer that is 
computationally and theoretically as simple as the dimer 
model, without introducing any approximation. 

In this paper we model the system in terms of a simple 
Hubbard, t—U, Hamiltonian, but the proposed procedure 
can be easily extended tot — U — Vort — J Hamilto- 
nians, to investigate charge-ordering transitions, whose 
possible coexistence with spin-order has been recently 
suggested.c3 Moreover, the symmetry properties can be 
conveniently implemented in more refined calculation 
schemes, to get simpler and more reliable description of 
the physics of K-phase salts. The paper is organized as 
follows. The next section is devoted to the description of 
the method. We then analyze the magnetic instabilities 
of ET-Y salts, and discuss the effects of the instabilities 
on the band structure. The difference between our sym- 
metrized mf and the other mf approaches is stressed, and 
the reliability of the dimer model is shortly addressed. 
Finally, we make connection with the experiment by dis- 
cussing the pressure dependence of the SC/AF interface 
in the ET-Cl, and by making a comparison with ambient 
pressure ET-Br superconductor. 
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II. THE SYMMETRIZED MEAN FIELD 
APPROACH 

Consistently witJa-experiniental data on centrosymmet- 
ric K-phase salts, E3IIJ we consider a unit cell with 4 equiv- 
alent molecular sites, and do not allow for modification 
of the periodicity of the crystal structure at the magnetic 
phase transition. We adopt the t — U Hubbard Hamil- 
tonian to describe Coulomb interactions giving rise to 
magnetic ordering: 

H= ^ Uj {al^aj^ +h.c.) + ^'^nini - U'^SiSi 

<i,j>(7 i i 

(1) 

where the indices run on the BEDT-TTF sites, the first 
term accounts for the intersite hopping, and the other 
terms describe the on-site Coulomb repulsion. In eq.(|^), 
Ui is the usual site number operator rii — ri^i -I- riii, and 
Si = ('T-iT ~ ^ii.)/2 is the net magnetization operator. 

In the mf approximation the many-body interaction 
is described by an effective single particle interaction, 
where each particle feels the other particles as a source 
of a mf potential. Then each product of two electronic 
operators AB is approximated with an expression where 
only a single operator appears, the effect of the second 
operator being substituted by its ground-state expecta- 
tion value. This approach gives reliable results when the 
fluctuations of the observables are small, although not 
zero as in single-particle approaches. Mathematically: 

Ab ^ {{A) +AA){{B) +AB) 

^ {A){B) + {A)AB + (B)AA (2) 

Thus in mf the two-particle Hubbard terms of eq. (|l|) 
become: 

i i i i 

(3) 

The equivalence of the four molecular sites imposes the 
constraint: 

+ K^} = 1.5 i = l..A. (4) 

Therefore the first term on the r.h.s. of eq. is a con- 
stant, and the relevant physics is described by the net 
magnetization term. We rewrite it by exploiting symme- 
try, and define within each unit cell the following four 
order parameters: 

VAFl = (si + S2 - S3 - S4) 

?7AF2 = (Sl - S2 - S3 + S4) 

VAF3 = (Sl - S2 + S3 - S4) (5) 

VFM = (si + S2 + S3 -I- S4) 



or, in short: r]^ = with v = AFl, AF2, AF3, 

FM. In these equations i counts the four BEDT-TTF 
sites within the unit cell, as indicated in fig. ^. The 
order parameters tjafi, Vaf2 and tjafs correspond to 
the three possible antiferromagnetic orderings; rjpM de- 
scribes the ferromagnetic phase. The net magnetiza- 
tion term in Eq. (|3|) then becomes: —2U^^{si)si = 

~{U /2)^- ^{ri^)riu^ where j runs over the unit cells. 
Since translational symmetry is not broken by the mag- 
netic transition, (77,^) is independent on j. The four mag- 
netic phases have different symmetry, so that the four 
•qi, order parameters are orthonormal and can be investi- 
gated separately, leading to the symmetrized mf Hamil- 
tonians: 

<l,k>,a j 

where tik are the hopping parameters, i.e., ti,^, ti,^, tp, 
and tq defined in fig. |^. For each symmetry, the ef- 
fective single particle potential Y^, is related to the ex- 
pectation value of the relevant order parameter by the 
self-consistency equation: 

Y.^^M (7) 




FIG. 2. Schematic view of the ac plane of centrosymmetric 
BEDT-TTF salts. 



By rewriting the last term in eq. (|6|) in terms of the orig- 
inal riii, riii operators, one immediately recognizes that 
Hi, is the sum of two independent tight-binding Hamilto- 
nians, H^i, describing electrons with up and down 
spin, respectively. The off-diagonal part of each one of 
these Hamiltonians is exactly the same as in the origi- 
nal tight-binding model, but the mf treatment of on-site 
electron-electron interaction introduces a diagonal contri- 
bution. Specifically, the diagonal elements of H^i^ within 
each unit cell are: 
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{H.^)u = -^< (8) 

and {H^-^)ii — —{H^i)ii. The two tight-binding prob- 
lems described by and H^i are easily diagonalized 
for different Yd values on very large lattices. In our ap- 
proach imposing the self-consistency relation on U sim- 
ply implies the ratioing of and {r]u), at variance with 
the lengthy and memory consuming iteratiopn^teps re- 
quired by a multi-parameter mf calculation.clLl This is 
very important in keeping the numerical procedure un- 
der control and allows us to work with very large lattices, 
tipically up to 6x10^ sites. Such large lattices, one order 
of magnitude larger than the largest lattice in Ref. ^, are 
diagonalized with no effort on a Digital Alpha 255 work- 
station equipped with 64 MB RAM. As we will discuss 
below, working on large lattices is very important to get 
an accurate description of the early stages of the phase 
transition, and then to get reliable information on the 
nature of the transition itself. 

The diagonalization of -ffi/j^ immediately defines 
the band structures for up and down spins. In the case 
of the FM instability, all the in eq. ^ are equal to 
1, so that, apart from a rigid shift of the energies by 
-{+)Yfm/'2 = -(+)[/ (?yFM)/4 for up (down) spins, the 
eigenstates are exactly the same as in the non-interacting 
case. Therefore, the originally degenerate bands for up 
and down spins are split by Yp^ = U{r]FM)- The Fermi 
level is fixed by the conservation of the total number 
of electrons, leading to unbalanced up and down spin 
population. If, without loss of generality, we consider 
positive (tjfm), we end up with lower energies for up 
spins and then with a ferromagnetic state characterized 
by larger population of up than down spins. 

In the case of AF order, instead, finite Y^, deform the 
original bands of the non-interacting system, due to the 
appearance of relevant diagonal terms in the real space 
Hamiltonian (eq. (||)). In this case, the eigenvalues for 
up and down spins are exactly the same, and the bands 
for the two spins stay exactly degenerate as in the non- 
interacting case, but the distribution of the two spin 
species is different on the sites, with a larger number 
of up spins residing on sites with negative coefficients 
Hvaf) > 0). 

III. RESULTS 

Table || summarizes the t's obtained from tiift-available 
structural data of ET-Cl and ET-Br sahsMH AU t's 
have been obtained from extended Hiickel (EH) calcula- 
tions on the pairs of BEDT-TTF molecules correspond- 
ing to the four interactions depicted in fig. ^ Specifically, 
they are evaluated as half of the splitting of the HOMO 
energy in each pair. It is well known that the values of 
the hopping integrals show large differenccp, depending 
on the method adopted for their estimate.E3 Therefore, 
the i's estimated for each structure and the resulting Uc 



have not to be assigned too much confidence. However, 
comparing results obtained with the same procedure on 
different structures is certainly informative. We have 
adopted EH estimates of since they compare well with 
available ab initio results.ES 
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FIG. 3. ET-Cl at T = 127 K, ambient p: t's from the 
first row of Table |. The order parameter for FM instability 
(squares) and AF instability of type 1 (circles) vs the effective 
mf potential Y , in the left panel, and vs U in the right panel. 
Filled circles mark points introduced in the linear regression 
of the rj vs. Y to determine Uc (see text), the arrow marks 
the point where the conduction gap opens, and whose band 
structure is shown in fig. Q bottom panel. 

Fig. ^ reports the Y^ and U dependence of (rjAFi) and 
(vfm) order parameters, as obtained for the t's relevant 
to ET-Cl at 127 K (first row in Table |). (?7af2) and 
iVAFs) curves are not shown since the comsponding in- 
stabilities occur at U larger than ^ 1 eVjlJ and are not 
relevant to our discussion. Indeed, even the ferromag- 
netic instability occurs at U higher than that for AF 
instabilities, and, in this respect, it is irrelevant from the 
physical point of view. However, the different behavior 
of (rjAFi) and {tifm) in the right panel of fig. |^ deserves 
some comments. 

Based on the two standard stability conditions: 
dE/drj — and d'^E/drf > 0, with E representing the 
expectation value of the working Hamiltonian, it is easy 
to prove that stable states for our system cojxpspond to 
points with positive slope in the ri{U) curvesE3 The neg- 
ative slope region in the rjFM(U) curve (marked by a 
dotted line in the right panel of fig. ^ then corresponds 
to unstable states, i.e. states that cannot be reached by 
our physical system. Thus the ferromagnetic instability 
corresponds to a first order phase transition, character- 
ized by a discontinuous jump of the order parameter at 
the transition, located at Uc ~ 0.945 eV. The small re- 
gion around Uc where the riFM{U) is non-single- valued 
corresponds to the hysteresis region, where two stable 
states coexist. 

The behavior of {rjAFij is different, with riAFi{U) hav- 
ing infinite slope at Uc ^ 0.64 eV. The infinite slope 
is a direct consequence of a strictly linear riAFiiXAFi) 
dependence in a fairly large region around the origin. 
In fig. |[ left panel, the filled circles show the points 
that fall on a single straight line, {rjAFi) = X^afi, 
with a squared correlation coefficient larger than 0.99998. 
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By applying the self-consistency condition in eq. (|^), 
one immediately gets steeply increasing tiafi values at 
a fixed U = Uc = 2/x ^ 0.639 eV. The coefficient 
X = d'^E/dY^p^ represents the susceptibility of elec- 
tronic system to the Yafi perturbation: the critical U 
is thus related to the inverse of the electronic suscepti- 
bility. 

Extracting Uc from the slope of the {t/afi) vs Yafi 
curve is a much safer procedure than searching for the 
minimum U where finite (tjafi) appears. The calcu- 
lated {rjAFi) values are affected by finite uncertainties, 
with a minimum intrinsic uncertainty given by the in- 
verse of the number of unit cells. Since (t?a_fi) enters the 
Hamiltonian matrix as a multiplicative factor for U , the 
uncertainty in {tiafi) implies an uncertainty in U , with 
6U/U = 6{riAFi)/ iVAFi)- Therefore, at small {tjafi) the 
relative uncertainty on U can be very large. This is by 
no means accidental, but reflects the intrinsic limitation 
of investigating phase transitions through finite size cal- 
culations. At the transition in fact the correlation length 
of the fluctuations in the order parameters are in princi- 
ple inflnite, so that calculations on finite lattices lead to 
large errors. 




FIG. 4. Band structure for ETCl, same parameters as in 
fig. ^. The upper panel refers to the non-interacting system 
(or equivalently to the system before the transition), the bot- 
tom panel corresponds to the point marked by an arrow in 
figM Y — 0.225 eV. The dotted fine marks the Fermi energy. 



energy difference between two extreme points in the two 
conduction bands (specifically between M-point in the 
upper band and X-point in the lowest conduction band) 
to measure the conductivity gap, A. Negative A implies 
overlapping bands and then metallic behavior, positive 
A measures the semiconducting gap. 

In summary, for the fs in the first row of Table |, 
relevant to ET-Cl at ambient pressure and T — 127 K, 
we observe a fairly sharp transition, at [4=0.639 eV, 
from a paramagnetic metal to an antiferromagnetic insu- 
lator, as shown by the semiconducting gap that opens 
up right at the transition (fig. ||). In our approach 
the metallic phase includes the superconducting state, 
since our Hamiltonian does not account for SC ccuipiling. 
The criticaUJ is similar to available experimcntalErca and 
theoreticaM estimates of the effective U in BEDT-TTF 
sahs, U ^ 0.5 - 1.0 eV. Therefore ET-Cl is just located 
at the metal/AF interface, in agreement with several ex- 
perimental observations (see below). Again, in agreement 
with experinp^iLt-and also with predictions of previous mf 
calculations ,ErBEj the AF phase is characterized by par- 
allel spins residing on the 1-2 dimer (fig.|^), as could also 
been inferred from simple arguments based on the dimer 
pictured 

Having developed a simple and efficient method to 
solve the mf problem for the ET-Y family, we can now 
play around with parameters trying to gain some infor- 
mation about the rich phase diagram of these systems. 
In fig. ^ the continuous lines show the [/-dependence of 
{VAFi) and A, calculated for the available t's relevant to 
ET-Cl at p = 3 and 27 kbar (Table |). The critical U 
increases with p from 0.64 eV at T = 127 K, ambient 
p, to ^ 0.68 eV and ~ 0.91 eV, at ambient T and p — 3, 
27 kbar, respectively. The increase of Uc corresponds to 
a stabilization of the metallic phase, and can justify the 
appearance of SC in ET-Cl under pressure, as we shall 
discuss in more detail in next Section. 
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It is interesting to investigate the evolution of the elec- 
tronic bands along the AF transition. Fig. ^ reports the 
band-structure calculated for the non-interacting metal- 
lic system, and for a system located just where the tran- 
sition goes to completion, i.e. the point marked by an 
arrow on fig. ^. The two conduction bands, that are 
partly overlapped in the metallic system, are split apart 
in the AF phase, opening a gap and then leading to insu- 
lating behavior. Fig. ^ reports the [/-dependence of the 
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FIG. 5. The conductivity gap (see text) vs U, same param- 
eters as in figj^. 

In the scale of fig. ^ the curve relevant to p — 3 kbar 
shows a very narrow region with a negative slope. How- 
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ever, the width of this region is only 2-3 times the nu- 
merical uncertainty on U, so that we cannot make any 
strong statement about observing a discontinuous, first 
order transition. In any case, the coexistence region, 
i.e. the hysteresis region for this transition, if present, 
would be so small to be irrelevant for any practical pur- 
pose. The region of negative slope disappears at p = 27 
kbar, where the transition looks smoother, with possibly 
a finite positive slope. Once more the effect is tiny and 
hardly disentangled from numerical uncertainties. 




A(eV) 



1.2 1.8 

U(eV) 

FIG. 6. The AF order parameters and the conductivity 
gap vs U. Continuous lines refer to t- values in Table |, with 
p increasing from left to right. Dashed lines refer to nominal 
pressures p — —33 kbar, at left, and p = 95 kbar, at right (see 
text). 



To get some clearer feeling about the role of pressure, 
we have linearly extrapolated the i-estimates available 
at p = 3 and 27 kbar to higher pressures. In fig. |^ the 
rightmost dashed line shows the corresponding evolution 
of the order parameter for a nominal p — 95 kbar, where 
tq extrapolates to zero. The smoothing of the transition 
is now evident: in this case {t]afi){U) clearly has a well- 
defined positive slope. An interesting observation is that 
the conduction gap (lower panel) closes not at the very 
beginning of the transition, but only when the transition 
comes to completeness. This corresponds to the appear- 
ance of a region of stability for a phase with simultane- 
ously AF distortion and residual metallic behavior. In 
other words, when the transition is continuous, our data 
suggest the presence of an antiferromagnetic metal phase. 



similar to that discussed at length in Refs. How- 
ever, the extent of this phase is very narrow, so that it 
represents at most a transient phase. 

Searching for some evidence of discontinuous phase 
transition, we have also "released the pressure" , by ex- 
trapolating the i's estimates at p = 3 and 27 kbar down to 
a nominal pressure of —33 kbar. Surprisingly, the transi- 
tion is continuous again, and smoother than at ambient 
p. Therefore, we find no clear evidence for discontin- 
uous transitions, except possibly in a very narrow re- 
gion in the parameter space corresponding to ET-Cl at 
p = 3 kbar. In general, the observation of very narrow 
regions of discontinuity and/or AF-metallic coexistence 
is strongly affectejd by numerical uncertainties and/or 
finite-size effectsr^ These transient states, being intrinsi- 
cally unstable, cannot be associated with physically sig- 
nificant states. On the other hand, their presence for 
particular values of the parameters signals an intrinsic 
instability of the system to external perturbations, such 
as those eventually leading to superconductivity. 



IV. DISCUSSION 

Overall our results agree witk-ntevious mf calculations 
on K-(BEDT-TTF)2Cu(NCS)2,fflt3 describing the tran- 
sition from a paramagnetic metal to an AF insulator oc- 
curring at Uc ~ 0.6-0.8 eV (the precise value of course 
depends on the choice of the t's). However, some de- 
tails on the evolution of the order parameter and on the 
opening of the semiconducting gap are different. In the 
recursive approach to the solution of the mf problem, 
adopted so far in the literature, the only viable proce- 
dure to estimate Uc relies on searching for the minimum 
U where finite {r]AFi) appears. As discussed above, this 
procedure leads to large uncertainties in Uc, that have 
to be properly accounted for ift-the analysis of numerical 
results. Kino and FukuyamauO use very small lattices 
(N=3600), corresponding to an intrinsic uncertainty in 
ivAFi) of at least 4 x 10"'*. In Ref. || the onset of AF 
is estimated to occur at Uc2, with {t]afi) ^ 0.02. This 
small value for the order parameter implies a minimum 
uncertainty in Uc2 of ~ 0.02 eV. Then the two transition 
points observed by Kino and Fukuyama, Uci =0.762 eV 
and Uc2 = 0.758 eV, coincide within numerical accuracy. 
In the lack of additional information their data are con- 
sistent with a single transition, as we find for the t's rele- 
vant to ET-Cl at T = 127 K or at p = 3 kbar (fig. ||) The 
presence of an antiferromagnetic metallic phase is then 
questionable. Similar problems occur in the interpreta- 
tion of data in Ref. |l^. Here a two-transition scenario is 
proposed at low p, involving a continuous transition from 
a paramagnetic metal to an antiferromagnetic metal (fi- 
nite rjAFi and negative A), immediately followed by a 
first order transition to an insulating state. This com- 
plex scenario, that we were unable to reproduce in our 
large lattice for any choice of the parameter set, is proba- 
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bly either a finite-size effect or a numerical artifact. One 
must also keep in mind the possibility it represents a 
characteristic feature of k(BEDT-TTF)2Cu(NCS)2, due 
to its lower symmetry. However, since it does not appear 
in the more symmetric ET-Cl phase at any pressure, it 
is irrelevant as far as SC is concerned. 

In Ref. H the numerical uncertainty in fixed 
by the authors at 0.001, propagates to give S{i]afi) ~ 
0.0014. As a consequence, the estimate for Ua, i.e. 
the critical U for the appearance of AF order, obtained 
for (rjAFi) = 0.01, is affected by a large uncertainty: 
Ua = 0.7 ± O.leV. More precise estimates are obtained 
for larger (rjAFi), e.g. Uc — 0.685 is essentially constant 
for {ijAFi) = 0.037, 0.499, 0.582, representing a good es- 
timate for the critical U where AF order appears and, at 
the same time, the electronic orbits close. The proposed 
estimate of the critical U for the opening of the semi- 
conducting gap, Ui = 0.699 ± 0.001 is different from Uc, 
again suggesting the presence of an intermediate phase 
with both metallic and AF character. Quite in agree- 
ment with our results at large p, the metallic antiferro- 
magnetic phase is a marginal phase that only survives in 
a very narrow transient regime. Indeed, as pointed out 
in Ref. |30| (Section 3.4.5), such a phase would imply a 
weak AF order and reconstruction of the Fermi surface, 
which however have not been experimentally observed. 
Our approach thus proves useful in excluding on a purely 
theoretical basis the spurious complexities in the phase 
diagram due to finite-size effects and/or numerical uncer- 
tainties, and should be particularly convenient when ex- 
tended to describe charge ordering instabilities together 
with magnetic instabilities. 

Several papers discuss K-phase BEEtLTTF salts 
within the the dimer approximation.QoE3o Basically, 
the tight-binding Hamiltonian for the four frontier molec- 
ular orbitals in the unit cell is rewritten in terms of the 
bonding and antibonding orbitals of the if,i-dimers. Since 
thi is at least twice as large as the other hopping inte- 
grals, the interactions between bonding and antibonding 
orbitals are neglected, and the original four-bands prob- 
lem reduces to a two-band problem. In the resulting 
lattice each site has four nearest-neighbor sites, interact- 
ing through [tp + tq)/2, and two next nearest-neighbors, 
interacting with t{,2/2. For the parameters relevant to k- 
phase salts, the bands calculated within the dimer model 
compare favorably with those obtained in the four-band 
calculation, confirming the validity of the dimer model 
approximation at least for the non-interacting case.E3 
The dimer-model lattice is simple, but still shows in- 
teresting physics. In fact, by varying the {tp + tq)/th2 
ratio, it interpolateSj-hetween a square lattice and a col- 
lection of ID chains.Eil Whereas it is suggestive to relate 
the variegated behavior obtained from such a model to 
the variety of observed properties for K-phase salts, some 
caution is in order. Just as an example, consider the case 
tb2 ~> 0, where the dimer lattice reduces to a half- filled 
square lattice with perfect nesting. As it is well-known, 
the critical U for the antiferromagnetic instability goes to 



zero in this limit, as also confirmed by mf calculations .113 
Instead, a mf calculation for the same parameters as in 
fig. 3, but ti,2 — 0, yields a continuous transition to the 
AF phase with a finite and fairly large Uc ~ 0.57 eV. This 
qualitatively different behavior can be easily rationalized: 
The small interactions between bonding and antibonding 
orbitals are large enough to break the ccuamensurability 
of the simple dimer model at tb2 = OEj'EJ A word of 
caution is also necessary when introducing electron cor- 
relations in the dimer model. Indeed, starting from a 
Hubbard Hamiltonian for the four-molecules layer, the 
resulting effective Udim for the dimer model is related to 
both U and tbi, according to a relation first proposed 
in Ref. ^ and rediscussed and extended in Ref. For 
the commonly accepted value of C/ ~ 4tfci, Udim is of 
the order of ttiB Thus the -applicability of FLEX and 
perturbartive approacheal3~tZl to the dimer model be- 
comes questionable, since these approaches work well in 
the limit Udim ~ ^bi ^ (^p + ^g)/2, ^62/2, where the dimer 
model itself breaks down. 

We now relate our results to the experimental observa- 
tions relevant to the ET-Y family. As mentioned above, 
several evidences indicates that the ET-Y salts are just 
at the AF/SC borderline. At ambient pressure the ET-Cl 
salt is a Mott antiferropiagnet, with a magnetic moment 
amplitude of 0.45 /iBj2i,thai compares favorably with 
the present and previougj^BO estimates of the magnetic 
order parameter. By applying pressure above 300 bar, 
ET-Cl shows a transition to complete superconductivity 
at about 12 K. At lower pressures, reentrant and partial 
superconductivity, with residual sample resistance, have 
been observed.E3 The fully deuterated ET-Cl (ds-ET- 
Cl) exhibits analogous behavior, only requiring a slightly 
higher pressure (440 bar) to reach superconductivity.L3 
ET-Br is superconducting at ambient pressure (Tc 11 
K) but the attainment of the superconducting phase is 
affected both by tiie cooling rate and by the deutera- 
tion of the sample.E3 It has been shown that by keeping 
constant the cooling rate at a sufhciently low value, the 
partially deuterated ^2- and (i4-ET-Br salts are super- 
conducting at practically the same Tc as the undeuter- 
ated sample. The de- salt, on the other hand, exhibits a 
complicated behavior attributed to the competition be- 
tween superconducting and insulating phase.L^ Finally, 
the fully deuterated dg-ET-Br is anj-antiferromagnetic 
insulator (magnetic moment: 0.3 hb)!^ and under pres- 
sure has a behavior similar to ET-Cl, jceaching complete 
superconductivity just above 60 bars.LJ For the sake of 
completeness, we mention that the ET-I salt is not su- 
perconclucting, even when pressures up to 5 Kbars are 
applicdB3 This kind of behavior has been ascribed to 
disorder. However, this compound is the least investi- 
gated in the ET-Y family, and we shall not consider it 
here. 

Rather obviously, a mf approach is inadequate to de- 
scribe phenomena related to non-equilibrium states, dis- 
order and/or sample inhomogeneity, like cooling rate ef- 
fects and reentrant superconductivity. We therefore focus 
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here on the complete-SC/AF crossover affected by pres- 
sure and/or by isotopic substitution. The universal phase 
diagram in fig. |^ presents the parameter "pressure" as the 
abscissa. It has been used to explain the differences in- 
duced by the Cl-Br substitution or dcutcration in ET-Y 
salts, and similar effects. In particular, the smaller radius 
of Cl~ Wth respect to Br^ implies a reduced effective 
pressuretJ in ET-Cl with respect to ET-Br (see Ref. ^ for 
a tentative numerical assessment of this effect). A similar 
effect of reduced pressure can be associated to deutera- 
tion, which corresponds to smaller end-group excursions 
around their equilibrium values. From the values of the 
corresponding Tc's, we can empirically associate an in- 
crease in p of ~ 380 bar for the Cl-Br substitution, and a 
decrease of ~ 140 bar for deuteration. A rationalization 
of these tiny effects is fairly difficult. We shall examine 
below whether our mf results help in this respect. 
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FIG. 7. Pressure dependence of the ratio between the areas 
of a and /3 orbits in the ac plane Fermi surface. The squares 
refer to the experimental data (Ref. ^5|), with the dashed line 
representing the best linear fit. Circles are the calculated 
values. 



We first focus on ET-Cl. The hopping integrals in the 
ac plane, calculated fos, the known crystal structures at 
arnbient p and 127 kJ23 and at 3 and 28 kbar (ambient 
T)t3 are reported in Table |. To make a first comparison 
with experiment, we have evaluated the areas of the a 
and /? orbits from the ac plane Fermi surface calculated 
in the tight binding approximation. The results relevant 
to the metallic phase are compared in fig. ^with the mea- 
sured areas from Shubnikov-de Haas (SdH) experiments 
performed at several pressureSria the 2-10 kbar range at 
the liquid Helium temperatureHj In order to renormalize 
the effects associated to the global volume contraction 
due to the different temperatures of the SdH and struc- 
tural data, the a//? ratio has been reported as a function 



of pressure. In fact, the /? area is equal to the area of the 
Brillouin zone for the 3/4-filled system. The 3 kbar aj (5 
nicely fits the experiment, and a straight line through the 
experimental points extrapolates near to the calculated 
28 kbar point. We notice that the crystallographic axes 
ratio c/a is practically unchanged with pressure.Ej This 
observation sheds doubts on the possibility of adopting 
c/a as a rough estimate of a//3, as suggested by Ref. ||. 

Table |, rightmost column, reports the ET-Cl Vc val- 
ues calculated by our mf approach. A word of warning 
is necessary when comparing data obtaiped at different 
temperatures, since it has been observedEa that the val- 
ues of the hopping parameters change correspondingly. 
Howe v£iV| this effect is not very pronounced in ET-Y 
family,Hc3 and we shall neglect it in the following. We 
notice that Vc increases monotonously with p, thus ac- 
counting for the pressure driven superconductivity tran- 
sition in terms of an increase of the critical value needed 
to reach the AF phase. In this respect, C/c seems to be a 
good "indicator" of the effective pressure of fig. 0. Ot 
previously suggested indicators, like the tbi/tp ratio| 
or the c/a ratio,0 seem to work less satisfactorily in this 
case: tbi/tp does not increase monotonously with p (Ta- 
ble^, and, as noted above, c/a is practically constant. 

We now turn attention to ET-Br. When comparing 
the ambient pressure, 127 K hopping integrals of ET- 
Cl with the corresponding ones of ET-Br, one finds small 
differences, and we get for the two systems virtually iden- 
tical phase transitions, occurring at basically the same Uc 
(Table |). Therefore, the different ground state of ET-Cl 
(AF) and ET-Br (SC) at ambient pressure cannot be un- 
derstood in terms of a difference in Uc- One might think 
that the actual effective U is different in the two types of 
salts, being smaller in ET-Br due a larger screening of the 
intersite Coulomb potential from the more polarizablc Br 
anions. However, this kind of qualitative explanation is 
not corroborated by the numerical values of the aajon po- 
larizabilities obtained from ab initio calculationsJ2ZI and is 
difficult to reconcile with the observation of an AF state 
for dg-ET-Br at ambient pressure. We could not calculate 
the hopping integrals in this the atomic coordi- 

nates are not available in the literature. We have used 
the f s calculated in Ref. |6|for both ET-Br and dg-ET-Br 
at 127 K (properly rescaled since the method of calula- 
tion is different from ours). We do not find significant 
difference between the UcS of the two compounds. 

As we already pointed out, c/a or ibi/tp are not good 
indicators of the properties of K-phase salts, and cannot 
be chosen as the x-axis parameter in a universal phase 
diagram like that reported in fig. |l|. Both the ratio of a 
and (3 orbits and Uc work satisfactorily as far as the p- 
dependence of ET-Cl properties is concerned. However, 
both fail if applied to rationalize the different behavior 
of ET-Cl and ET-Br and/or the effects due to deuter- 
ation. It is important to underline that both a/ (3 and 
Uc are "single particle" parameters, in the sense that 
they are fully determined by the band structure, i.e., the 
t's values. Investigating the band structure of K-phase 
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salts offers useful information to rationalize their behav- 
ior, but this information is not enough, and the role of 
interactions beyond single-particle picture has to be in- 
voked to understand the behavior of systems near the 
AF/SC interface. The failure of simple band-structure 
treatments for K-phase salts has been recently suggested 
based on the padependence of cyclotron effective masses 
with pressurej£3 and can also be recognized from high re&a 
olution mesurements of thermal expansion coefficients.ES 
Residual electronic correlations, disorder induced local- 
ization effects, electron-phonon coupling and interlayer 
effects all can play an important role, particularly at 
the AF/SC interface. More theoretical and experimental 
work is in order to settle the relative importance of these 
effects. 
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TABLE I. Hopping integrals and critical U for ET-Cl and ET-Br salts. All parameters in eV. 

^61 ib2 tp tq Uc 

ET-Cl, amb. p, T =127 K 0.2315 0.0760 0.0901 0.0410 0.639± 0.001 

ET-Cl, p=3 Kbar, amb. T 0.2239 0.0851 0.0844 0.0517 0.676± 0.005 

ET-Cl, p=27Kbar, amb. T 0.2770 0.0935 0.1400 0.0380 0.906± 0.005 

ET-Br, amb. p, T =127 K 0.2244 0.0712 0.0936 0.0396 0.636± 0.001 
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